Path to finding tlie critical thickness for memory in thin ferroelectric films 
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The finite screening length by real metallic electrodes, albeit very small (< lA), results in finite 
depolarizing field that tends to split the film into domains. In very thin ferroelectric films the 
domain structure reduces to sinusoidal distribution of polarization considered first in the 1980s. We 
discuss the phase transition between this structure and a single domain state and show that it is 
first order, if it exists at all. The alternative possibility is that the single domain state at zero bias 
voltage would be metastable for all temperatures in most systems. This scenario defines a path 
towards solution to a problem of finding parameters of a system that can sustain the ferroelectric 
memory over a desired period of time. 

PACS numbers: 77.80.Dj, 77.55. +f, 77.22.Ej 
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Stability of ferroelectricity in ultrathin films with elec- 
trodes is a topic of intense current interest, as illustrated 
by numerous papers that appeared recently in the leading 
journals (see, e.g. flSiHi]). In this Letter, we show 
that all those papers have little or no relevance to the 
possibility of the ferroelectric (FE) memory since they 
do not actually address the core issues. We shall formu- 
late the actual problem that has to be resolved in order to 
answer the question at hand, and indicate possible paths 
to its solution. Our present formulation builds on an old 
paper by Chensky and Tarasenko (ChT) [s], which, un- 
fortunately, is rarely mentioned in recent literature. ChT 
discussed stability of the FE films with respect to small 
fluctuations of the single and multidomain states in an 
electroded film of a uniaxial ferroelectric with 'dead' lay- 
ers between the electrodes and the material. Fig. la (in- 
set). It was already shown that their system satisfactorily 
models also the FE films without a physical dead layer 
but real metallic electrodes Q, i-e. it considers the same 
type of a system as have the recent Refs.P, S S IHIj 
which customarily disregard the multidomain region. 

The results of ChT and the above authors for stability 
of homogeneous FE state in the films would only be cor- 
rect if the phase transitions between the single and mul- 
tidomain states were second order. In fact, the character 
of the transition between the ferroelectric single domain 
state and the multidomain ones was not studied by ChT, 
and it is one of the main questions we address in this Let- 
ter. We find that in the low-T part of the phase diagram 
the transition is first order. Thus, the stability of the FE 
state found by ChT means, in fact, its metastability at 
least in some temperature region. With account of our 
old result Q, we expect that the single domain state is 
metastable for all temperatures or, if one studies the FE 
states at a given temperature but for different film thick- 
nesses, at all film thicknesses. This question has to be 
studied in more detail, but in any case, the key theoreti- 
cal problem to be resolved in order to predict feasibility 
of the FE memories is not that of existence of a solu- 



tion corresponding to the ferroelectric phase or demon- 
strating its stability with respect to small fluctuations. 
For the case of metastability of a single domain state, 
the key problem is calculating the escape rate from this 
metastable state. We are not aware of theoretical efforts 
to solve this problem. Another aspect of the memory 
problem is to find the conditions of absolute stability of 
the single domain state under zero external bias voltage. 



The single domain FE state is evidently stable in the 
case of the ideal metallic electrodes (here and below we 
mean, of course, the ideally homogeneous systems.) For 
a non-zero dead layer thickness {d) this is not evident 
anymore, although it is hardly surprising that there ex- 
ists a "minimal thickness" dm such that for d < d„i the 
FE single domain phase is stable in the ChT sense, i.e. 
with respect to small fluctuations, in all the temperature 
range. The minimal thickness dm depends on parameters 
on the FE material, as well as on the dielectric constant 
of the dead layer (or, in the case of real metallic electrode, 
on parameters of both the electrode and the ferroelectric 
film.) The stability in the ChT sense of single domain FE 
state does not mean that this state is absolutely stable. 
According to our result in Ref. [Sj, at any finite thick- 
ness the homogeneously polarized state is metastable far 
enough from the phase transition (or for thick enough 
films.) However, it is stable not very far from the tran- 
sition or for very thin films. This conclusion makes this 
case attractive from the point of view of memory applica- 
tions. Unfortunately, the experimental systems of today 
do not fall into this category It does not seem im- 
possible, at least in principle, to find a suitable system, 
but we shall not dwell on this issue here, assuming below 
that d > dm- 



Similar to ChTjej, we consider a uniaxial ferroelectric 
film under the external bias voltage U (Fig. la, inset). 
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The free energy of such a film is given byji 
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where Fq is the free energy of the system at P = 0, 
with Pz the ferroelectric (switchable) component of po- 
larization, Pzb the nonferroelectric part of the polariza- 
tion perpendicular to the electrodes Q, P± the in-plane 
polarization, A = j {T — Tc) , and B,Dij,r] = const, 
Vj^= {d^,dy) is the gradient in the plane of the film, 
i,j = x,y, and we assume summation over repeating in- 
dices, A± > 0, E the electric field, Ce the dielectric con- 
stant of the electrode/dead layer (marked DL), Q w —p 
the electrode charge, p = (Pz) is the average polarization 
in the FE film. We assume that Dij = DSij, which is 
valid, in particular, for BaTiOs or PbTiOa films grown 
on (100) SrTiOs substrate because of a square symmetry 
in the film plane. The noncritical in-plane and out-of- 
plane dielectric constants are equal e± = I + AttP^/ E± = 
1 + ATr/A± and €{, = 1 + Att/k, respectively. 

The stability of the paraelectric phase is lost with re- 
spect to appearance of the "polarization waves" @, Q : 



Pz (x, z) = a cos qz cos kx, 



(2) 



where, for d > 2d„ , g ~ ir/l and k = {Att^ /ej_DpY^*(6, 
0] . To find the amplitude of the polarization wave in the 
external bias field, we use the following approximation 
valid close to the phase transition (see below): 



Pz{x, z) = p + P = p + a cos qz cos kx. 



(3) 



In this case, the non-equilibrium free energy per unit area 
F{p, a) takes the form: 
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where A = A + 2DP = 7 (T - r<j) , ^ = 
And/ (eJ + Cbd) ~ 2Dk^ « And/ {ej) - 2Dk^, Eq = 
Cell / {ej + efjd) K, U /I the external field for the usual case 
of a thin dead layer eJ ^ ebd- At ^ > (or d > dm)^ 
according to the above expression for the free energy, 
the system will undergo a phase transition at Pq = 
and A = j/ = 0,i.e. aX T = Td into the sinusoidal 
domain phase. At lower temperatures, there may be 
another transition into a homogeneously polarized state 
with a = 0, p 7^ 0. The above potential allows one to 
study both phase transitions also at non-zero external 
field. 

When writing Eq. ([3]), we have assumed the stripe- 
like sinusoidal domain structure. This is far from being 



obvious in our isotropic case, and ChT discussed pos- 
sibilities of checkerboard and hexagonal sinusoidal do- 
main structures. However, those become irrelevant if 
one takes into account the elastic strains. Their cou- 
pling to the inhomogeneous polarization produces anhar- 
monic terms in Eq. @ , with the renormalized coefficients 
that now depend on the direction of the "wave vector" k 
for any elastically anisotropic medium [l3|- For the case 
of tetragonal uniaxial ferroelectric, the square symmetry 
tells us that there are at least two orthogonal orienta- 
tions of stable stripe structures with the same free energy, 
while a hexagonal domain pattern is clearly impossible. 
The checkerboard domain lattice should be studied sep- 
arately, but it is unlikely to be relevant. In the following, 
we discuss a stripe structure neglecting the elasticity and 
assuming that its role reduces mainly to selecting the 
direction of the sinusoidal polarization waves. To be pre- 
cise, elastic coupling also leads to changes of the coeffi- 
cient i? in (III), which are slightly different between the 
third, fourth, and fifth terms there (see Ref.[l3l), but this 
is only a numerical difference that does not affect any of 
the results below, and we shall not dwell on this issue 
here. 

By minimizing the free energy ^ with respect to a, 
one finds the equilibrium amplitude of the wave: 
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where s = p/pc, with pc — \J —A/iB the characteristic 
polarization. We finally arrive at the dimensionless free 
energy, / = iB^l-^F: 
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iy^s* - Vyse, \s\ < 1, 



(6) 

where e = Eq/C the relative external field, C = ^'^/^/VSP 
the characteristic electric field, and 



y = -A/i, 



(7) 



the characteristic temperature (i.e. the relative distance 
of transition temperature from the paraelectric phase 
that depends on the film thickness /). It is easy to see 
that the free energy is continuous with the first derivative 
with respect to s, while d^/_/ds^ ^ d^/+/c?s^ at s = ±1. 

One of the easiest ways to reveal the order of the tran- 
sition is to inspect the equations of state curves s = s (e) 
obtained from the condition df/ds = 0, which read: 
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Isl < 1. 
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FIG. 1: The equation of state p = p{Eo) for ferroelectric film 
with dead layer in external field Eq for various values of the 
relative temperature y, Eq. (a) y — 1/4, (b) y — 1, and 
(c) y = 3. There is a phase transition between the homoge- 
neous state and the one with sinusoidal domains: (a) second 
order, which becomes first order in cases (b) and (c), where 
arrows indicate the hysteretic behavior during polarization 
switching. The terminal external fields for hysteresis are at 
points A and B. All those cases are very different from the 
FE film without the dead layer. ( = £,^^^/VW is the char- 
acteristic electric field. Inset (a) shows the schematic of the 
ferroelectric film with either real electrodes or the dead layer 
with the thickness A. 

Recall that the given state is (ineta)stable only when 
cPf/ds^ > 0. Several typical s (e) curves for y — 1/4, 
1, and 3 are displayed in Fig. 1. We see that at a rel- 
atively small y — 1/A (not far below Td or A = y = 0) 
and in the field e > Cc = y/y — 2?/'^/^/3 = 5/12 the sys- 
tem is homogeneously polarized. In the lower bias field, 
it splits via the second order phase transition into do- 
mains with zero net polarization at £^o = 0- At y > 3/8, 
the transition is first order, which is clearly seen for 
y = 1, where two metastable solutions exist in the fields 



CA < Eq/( < e_B, Fig. 1. The points of the first order 
phase transitions found from the condition of equal free 
energies, /+[s+(e)] = /_[s_(e)], are marked in Fig. 1. 
The first instance when the state with spontaneous net 
polarization p ^ ai Eq = becomes formally possible 
as a solution to the equations of state is at ?/ — 3/2. How- 
ever, this state is unstable {(Pf-/ds^ < at s = ±1), as 
is evident from the negative slope of the lower branch of 
p ~ e curve at s = 1 in Fig. lb in the external fields 
eA < -E-o/C < es- The above tricritical behavior of 
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FIG. 2: External field-temperature (thickness) phase dia- 
gram of a ferroelectric thin film with real metallic electrodes. 
The line of stability loss of the homogeneous states {E^, E^) 
has been calculated in Ref. [^. The tricritical point (y = 3/8) 
and lines Ei (line of first order phase transitions), Eb are 
from the present work. The point yi = 1 corresponds to the 
"critical thickness" calculated in Refs. P, S H 0] with an 
account for atomistic structure of the electrodes (additional 
boundary conditions), j/2 = 3/2 is the stability point consid- 
ered in Ref. with an account for elastic strains. The shaded 
region is area of interest for estimating the retention time of 
homogeneously polarized state. 

thin ferroelectric films with real electrodes becomes clear 
from the phase diagram in {Eo,y) plane. Fig. 2. Indeed, 
there is a second order phase transition in the interval 
< y < 3/8, where the sinusoidal domains form in the 
fields below e = ca, where 

eA = Ea/C = - 2y/3), (10) 

corresponding to point A in Fig. lb. In the range 
3/8 < y < 3/2, the same expression gives the line of sta- 
bility loss of the single domain phase. After passing the 
tricritical point at y = 3/8, there appears a line of first 
order phase transitions Ei(y), that terminates at y = 3, 
and it is bracketed by the lower and higher terminal fields 
CA (fTU|) and cb, where 

eB^EB/C = 2{3 + yf/y 27. (11) 

The physical meaning of these fields is clear from Fig. 2: 
Ea is the lowest field where the stability of the homo- 
geneous phase is lost with respect to domains, while Eb 
is the highest field where the domain structure still ex- 
ists. The hysteresis, therefore, is observed in the range of 
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fields Ea < Eq < Eb, when y > 3/8. Although the ap- 
proximation of sinusoidal domains becomes rather poor 
at larger parameters y ^ yi = 1, it should still correctly 
grasp the main features of the phase behavior. Note that 
this result invalidates a long unchallenged claim by Chen- 
sky and TarasenkoQ that one can prepare a monodomain 
state at E'o = 0, y > ?/2 = 3/2 (that corresponds to low 
temperatures well below Tc and Td) by polarizing the 
system in sufficiently high field and then removing the 
field. 

The present diagram suggests that the homogeneously 
polarized state will remain metastable (in the present 
one-sinusoid approximation) only at y > 3. There is a 
metastability of homogeneously polarized state in the re- 
gion 3/2 < y < 3. Formally, both the ferroelectric mem- 
ory and the polarization switching are possible at these 
temperatures/thicknesses but no conclusion of practical 
importance can be made before calculating the escape 
time from the metastable state. At larger y > 3 (further 
down from the phase transition with respect to temper- 
ature, or for films with thickness exceeding the critical 
one) the state with the homogeneous polarization in the 
present single- harmonic approximation, Eq. ([3]), has the 
same or lower free energy than the state with p — : 
f+{sm) ^ /-(O), where ±Sm are the positions of the min- 
ima of the free energy /+ at iJg = (Fig. 3). Note that 
this result is approximate. The reason is that Eq. ([3]) is 
valid near the phase transition point only. The region of 
validity of this approximation has been estimated in 
as roughly —A < Dk^, which means y < 1 if 47r(i/ (ej) 
and 2Dkl are of the same order of magnitude. 

We should mention that the second order phase transi- 
tion into homogeneous FE phase, considered in the prior 
papers [E S S S|j may only occur in zero field, Eq = 0, 
and it corresponds to A = And/ {ej + Cbd) , or, in other 
words, to the point y = yi = 1 on the {EQ,y) diagram, 
which is the point that in typical circumstances lies well 
inside the domain regime and, therefore, is normally 
unreachable, Fig. 2. On the other hand, Pertsev and 
KohlstedtQ have studied the stability loss of the FE ho- 
mogenous state, which corresponds to point y = y2 = 3/2 
in the phase diagram. They noted that one should take 
elastic coupling into account while discussing the stabil- 
ity of the ferroelectric state but, unfortunately, did it in 
a confusing manner with incorrect conclusions, see anal- 
ysis in In a more accurate approximation account- 



ing for higher harmonics to describe the inhomogeneous 
polarization, the free energy minimum at s = ii'o = 
dips lower than that of the homogeneous state. It is 
those higher harmonics that convert the sinusoidal do- 
main structure into a conventional one with narrow do- 
main walls. For the parts of the curves corresponding to 
|s| > 1, these higher harmonics are not important (they 
are when an amplitude of the first harmonic becomes sub- 
stantial) but they will change the curves for |s| < 1 sub- 
stantially. The amplitudes of the higher harmonics are to 
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FIG. 3: The free energy of the FE film with the dead layer for 
the relative temperature y = 3, when the transition between 
the homogeneously polarized state and sinusoidal domains is 
first order. It is evident that switching proceeds through the 
state with domains with much lower energy barrier for nucle- 
ation (top of the barriers indicated by arrows). An account for 
higher harmonics in Eq. ([2]) will deepen the energy of the do- 
main state compared with the one shown here and will lower 
the barrier even further. 



be considered as new variational parameters for the free 
energy, and their account will be lowering the estimated 
free energy. Hence, the minimum at s = p = in Fig. 3 is 
actually deeper, and the homogeneously polarized phase 
becomes stable not at y = 3 but at a larger value (i.e. at 
a lower temperature or a larger thickness.) Furthermore, 
it is possible that the homogeneous state would always 
remain less stable than the polydomain state in that re- 
gion. Indeed, in the opposite hmiting case, i.e. far below 
the FE transition, it has been shown that for any thick- 
ness of the dead layer the multidomain state has lower 
free energy than the homogeneously polarized state [1]. 
Anyway, what value of y would correspond to the "criti- 
cal thickness for the ferroelectricity" within the discussed 
simple case depends on the desirable memory retention 
time and should be found by solving a kinetic problem. 

The present discussion of an equilibrium problem in 
the one-sinusoid approximation can be improved. We 
have already mentioned the effect of coupling of the po- 
larization and the elastic strains, which is unimportant 
for studying the stability of the paraelectric phase but 
should be taken into account while considering the mul- 
tidomain state. This is valid also for the discussion of 
stability of a single domain FE phase. We can men- 
tion also the effect of the additional boundary conditions, 
apart from electrostatic ones used by ChT, and the ef- 
fects of higher order terms in the LGD expansion. This, 
however, would not change the present qualitative con- 
clusions that will apply also to finite FE patches with 
lateral sizes of lOOnm and smaller, which have lateral di- 
mensions still much larger than the domain width, which 
is just ~ 1 — 2nm in 5nm thick BaTi03 film0]. Our dis- 
cussion above indicates only the first steps on the path 
to addressing the kinetic problem in a likely scenario 
that the single domain state is not absolutely stable but 
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